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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted 
(metric)  units  as  follows; 

_ tfaUiyly _  _ Sy - 

cubic  feet  0.02831685 

degrees  (angle)  0.01745329 

feet  0 . 3048 


to  SI 

To  Obtain 
cubic  metres 
radians 
metres 
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DEPTH-AVERAGED  NUMERICAL  MODELING 
FOR  CURVED  CHANNELS 


PART  I :  INTRODUCTION 

Genuinely  two-dimensional  (2-D)  flows  are  uncommon  in  hydraulics,  but 
depth— averaged  2— D  models  have  been  widely  used  for  shallow  water.  One 
shortcoming  of  such  models  is  that  they  fare  poorly  in  channels  with  signifi¬ 
cant  curvature.  This  is  due  to  the  models*  inability  to  account  for  the 
effects  of  helical  secondary  flow,  which  is  a  three-dimensional  (3-D)  phenome¬ 
non  that  cannot  be  represented  by  conventional  depth  averaging. 

When  the  depth-averaged  streamlines  are  curved,  and  the  vertical  distri¬ 
bution  of  velocity  is  nonuniform,  centrifugal  forces  create  a  torque  that 
generates  hellcity  in  the  flow.  (This  is  like  a  screw  advancing  or  retreating 
in  the  streamwlse  direction. )  If  no  correction  is  added  for  the  secondary 
flow,  depth-averaged  models  will  usually  constrain  the  largest  velocities  to 
lie  near  the  inside  of  a  channel  bend,  instead  of  allowing  the  observed 
gradual  migration  of  high  velocity  toward  the  outside. 

The  influence  of  secondary  flow  is  well-known  for  natural  meandering 
channels,  where  it  alters  the  depth-averaged  flow  and  ultimately  the  bed 
itself.  In  the  past  this  has  been  simulated  with  both  analytical  (closed- 
form)  and  numerical  (discrete)  models.  Especially  noteworthy  among  the  ana¬ 
lytical  models  is  the  work  of  Johannesson  and  Parker  (1989a),  who  developed 
eqtiations  for  predicting  velocity  redistribution  in  meandering  rivers  (1989b). 
Their  success  arises  from  two  Important  model  features.  First,  a  perturbation 
expansion  was  used  to  derive  a  differential  equation  for  the  secondary  flow 
from  the  3-D  equations  of  motion  (with  an  empirical  shape  function  for  the 
vertical  distribution  of  primary  velocity) .  The  differential  equation  was 
accurate  only  to  first  order,  but  its  closed-form  solution  offered  a  credible 
representation  of  the  secondary  flow  for  small  curvature  (large  radius). 
Second,  by  accounting  for  the  rapid  attenuation  of  all  velocity  components 
near  channel  sidewalls,  the  model  reproduced  the  velocity  redistribution 
observed  in  channels  with  uniform  width  and  variable  curvature. 

Depth-averaged  2-D  models  such  as  the  STREMR  code  (Bernard  1989,  1992) 
are  often  employed  for  shallow  channels  with  irregular  lateral  boundaries  as 
well  as  nonuniform  depth  and  curvature.  This  kind  of  variability  makes  it 
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inpractical  to  seek  a  closed— form  solution  to  the  Johannesson-Parker  equation 
for  secondary  flow,  and  It  may  even  stretch  the  applicability  of  the  equation 
Itself.  For  2-D  models  that  seek  only  to  calculate  the  depth-averaged  primary 
flow,  however,  the  overall  effect  of  the  secondary  flow  is  more  important  than 
its  actual  stinicture.  In  practice  it  may  be  sufficient  to  solve  an  empirical 
transport  equation  for  streamwise  vorticlty  that  reproduces  the  depth-averaged 
influence  of  the  secondary  flow,  but  not  necessarily  the  underlying  details. 
This  assumption  is  the  basis  for  the  approach  taken  herein. 

From  qxialitatlve  arguments  concerning  the  interaction  of  lateral  curva¬ 
ture,  bottom  friction,  and  vertical  nonuniformity,  one  can  establish  a  plausi¬ 
ble  form  for  an  equation  that  governs  the  production,  transport,  and  dissipa¬ 
tion  of  streamwise  vorticity  (helical  secondary  flow) .  The  latter  creates  a 
shear  stress  in  the  vertical  plane,  parallel  to  the  primary  direction  of  flow. 
Sidewalls  and  variations  in  depth  cause  lateral  nonuniformity  of  this  stress, 
which  transports  momenttun  in  the  lateral  direction  and  alters  the  distribution 
of  depth-averaged  velocity.  Empirical  coefficients  allow  the  equations  to  be 
ttoied  for  q\iantitative  agreement  with  experimental  data.  If  the  model  is  to 
be  truly  useful  for  making  predictions,  then  the  empirical  coefficients  must 
be  universal ;  i . e . ,  one  set  of  values  should  be  adequate  for  all  channel 
configurations  of  interest. 

This  report  documents  the  empirical  development  and  implementation  of  a 
secondary  flow  correction  (SFC)  for  the  2-D  STREMR  code.  Part  II  discusses 
the  depth-averaged  equations  for  the  primary  flow  and  the  k-e  turbulence 
model,  and  Part  III  presents  arguments  used  in  the  formulation  of  the  SFC. 

Part  IV  offers  a  comparison  of  computed  and  measured  results  for  three  differ¬ 
ent  bendway  configurations,  and  Part  V  summarizes  the  conclusions  and 
recommendations  for  this  study. 
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PART  II:  PRIMARY  FLOW 


The  governing  equations  for  the  primary  2-D  flow  are  the  depth-averaged, 
incompressible  Navier-Stokes  equations.  These  are  the  equations  for  conserva¬ 
tion  of  isass  and  momentum,  given  respectively  by 

V*(hu)  =  0  (1) 


^+uxVu=T-X  +  S-  ZE 

It  -  -  -  -  -  p 


where 

V  -  gradient  operator 
li  -  velocity 
t  -  time 

I  -  viscous  force  arising  from  the  depth-averaged  stress  tensor 
X  -  bottom-friction  force 
S  -  force  arising  from  secondary  flow 
p  -  pressure 
fi  -  density 

In  the  above  equations,  an  underline  indicates  a  vector.  The  cartesian  x-  and 
y-components  of  1  are,  respectively 

Tj  =■  h‘^i/V«{hVu)  ♦  2j/jjU3j  +  i/y(Uy+v^) 

Tj  ■  h'^|/V*{hVv)  2|/yVy  +  |/3j(Vjj+Uy) 

where  u  and  v  are  the  x-  and  y-components  of  u  ,  and  the  subscripts  x 
and  y  indicate  spatial  derivatives.  The  quantity  i/  is  the  eddy  viscosity, 
which  is  a  kinematic  viscosity  that  is  used  to  represent  the  influence  of 
small-scale  turbulence.  It  is  related  to  the  turbulence  energy  k  and  the 
turbulence  dissipation  rate  e  by 


V  ■  C„ 


idiere  is  a  dimensionless  empirical  coefficient. 

The  streamwise  deceleration  X  due  to  bottom  friction  is  represented  in 
STREMR  by 


6 


X  =  C,h-iu|u|  (6) 

where  Is  the  friction  factor  given  by  Manning's  equation, 

C£  =  9.81n2h'^/3  (7) 

and  the  quantity  n  is  Manning's  coefficient.  The  value  of  9.81  for  the 
constant  of  proportionality  applies  only  when  h  is  given  in  metres.  If  h 
is  given  in  feet,  this  value  becomes  14.5  Instead. 

The  streamwise  acceleration  S,  arises  from  nonuniformity  in  the  depth- 
averaged  shear  stress  r,  created  by  the  secondary  flow.  This  is  given 
approximately  by 


S  »  p-i^[h-^n«V(hr,)  +  2r-ir,] 


where  s  is  a  unit  vector  normal  to  u  ,  and  r  is  the  lateral  radius  of 
curvature,  which  is  related  to  the  velocity  and  its  derivatives  by 


UV(Vy  -  U,)  +  -  V^Uy 

The  operator  D’V(hr,)  is  the  normal  derivative  of  the  depth-integrated  shear 
stress  created  by  the  secondary  flow. 


n»V{hT,)  = 


v(hr4  -  u(hr.) 


The  k-«  turbulence  model  (Launder  and  Spalding  1974)  provides  semi- 
empirical  governing  equations  for  k  and  c  , 

♦  u»Vk  »  i/r  -  «  +  a^-ih-W»(h»/Vk)  ( 


2 

+  u»Vc  =  c^vV L  -  C2^  *  <7j*^h'^7»(hj/V£) 


In  Eqxiations  11  and  12,  the  second  term  on  the  left  is  the  advection  term. 
The  second  and  third  terms  on  the  right  are  the  dissipation  and  diffusion 
terms,  respectively.  The  first  term  on  the  right  is  the  production  term, 
which  is  proportional  to 
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r  =  2(u,2+Vy2)  +  {Uy+V,)2 


(13) 


In  this  context,  advection  means  transport  by  the  primary  flow;  produc¬ 
tion  means  creation  of  small-scale  turbulence  from  the  primary  flow;  dissipa¬ 
tion  means  frictional  loss  through  the  smallest  eddies  of  the  turbulence;  and 
diffusion  means  the  spreading  that  occurs  because  of  eddy  viscosity.  The 
standard  set  of  dimensionless  empirical  coefficients  (Patel,  Rodi,  and 
Scheurer  1985)  is 

C„  -  0.09 
Cl  -  1.44 
C2  -  1.92 

Ok  -  1.0 

CT,  =  1.3 

With  suitable  boundary  conditions  for  u,  v,  p,  k,  and  e,  and  with  suit¬ 
able  approximations  for  X  and  S  ,  Equations  1,  2,  11,  and  12  are  sufficient 
for  calculating  depth-averaged  primary  flow  within  the  limitations  of  the  k-e 
turbulence  model  and  the  definitions  for  X  and  S  . 

In  the  STREMR  code,  velocity  components  normal  to  the  boundaries  are 
held  fixed  on  inlets  and  solid  walls,  and  computed  by  a  discrete  radiation 
condition  (Orlanski  1976)  at  the  outlets.  The  total  flow  rate  remains  con¬ 
stant,  as  do  the  individual  flow  rates  through  each  continuous  inlet  and  out¬ 
let.  In  a  given  time-step,  the  velocity  normal  to  any  boundary  segment  is 
either  constant  (for  inlets  and  solid  walls)  or  determined  by  neighboring 
velocities  in  the  previous  time-step  (for  outlets).  For  calculating  vorticity 
and  shear  stress  on  sidewalls,  tangential  velocity  is  assumed  proportional  to 
the  1/7  power  of  distance  from  the  walls. 

The  value  of  k  is  fixed  at  a  small  fraction  (0.006)  of  the  primary 
inflow  energy  at  inlets,  and  the  normal  derivative  of  k  is  set  to  zero  at 
outlets .  This  estimate  for  the  inflow  turbulence  energy  (whose  actual  value 
is  unknown)  may  be  somewhat  low,  but  the  error  has  little  effect  downstream 
because  k  grows  rapidly  with  the  developing  flow.  On  the  sidewalls,  k  is 
assumed  proportional  to  the  1/7  power  of  distance  from  the  walls.  The  normal 
derivative  of  e  is  set  to  zero  on  all  boundaries  except  inlets ,  where  the 
value  of  €  is  fixed  to  maintain  a  Reynolds  number  of  50  based  on  the  eddy 
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viscosity,  the  streamwlse  grid  spacing,  and  the  Inflow  velocity. 

STREMR  uses  a  staggered  marker-and-cell  grid,  with  normal  velocity  com¬ 
ponents  defined  on  the  cell  faces  and  pressures  defined  at  the  cell  centers. 
This  finite— volume  grid  arrangement  allows  the  pressure  to  be  computed  from  a 
Poisson  equation.  Since  the  normal  velocity  components  are  known  (or  calcu¬ 
lated)  in  advance  for  all  boundaries ,  the  derivative  of  pressure  normal  to  the 
boundaries  can  be  set  to  zero.  STREMR  uses  a  variant  of  the  MacCormack  pre¬ 
dictor-corrector  scheme  (MacCormack  1969;  Bernard  1989,  1992)  to  solve  Equa¬ 
tion  2,  and  a  single-step  upwind  scheme  (Anderson,  Tannehill,  and  Fletcher 
1984)  to  solve  Equations  11  and  12.  Steady-state  solutions  are  achieved  by 
marching  forward  In  time,  with  potential  flow  or  some  other  mass-conserving 
flow  used  for  the  initial  velocity,  and  small  uniform  values  specified  for  the 
initial  turbulence  quantities. 

For  open  channels,  the  water  surface  is  assumed  to  be  a  rigid  lid.  With 
this  constraint,  the  pressure  computed  for  any  given  cell  is  equal  to  the 
displacement  that  would  otherwise  occur  for  a  free  surface  at  that  location. 
This  assumption  limits  the  applicability  of  STREMR  to  subcritical  flow.  If 
the  streamwlse  variations  in  channel  width  or  depth  are  very  gentle,  and  if 
there  are  no  obstacles  (piers,  dikes,  or  islands)  in  the  flow,  then  the  code 
may  be  valid  at  Froude  numbers  as  high  as  0.6  or  0.7.  Otherwise,  the  upper 
limit  on  the  Froude  number  is  0.5. 
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PART  III;  SECONDARY  FLOW  CORRECTION 


To  obtain  an  approximation  for  r,  without  a  fully  3-D  discretization, 
one  first  needs  an  equation  for  production,  transport,  and  dissipation  of 
streamwlse  vorticity.  Qualitative  arguments  dictate  the  form  of  this  equa¬ 
tion,  whose  empirical  coefficients  (free  parameters)  have  to  be  adjusted  to 
make  model  predictions  fit  experimental  data.  If  the  SFC  is  to  be  quantita¬ 
tively  useful  in  general,  then  coefficients  obtained  for  one  channel  should 
yield  acceptable  predictions  for  other  channels  with  different  depth,  curva¬ 
ture,  and  bottom  friction. 

Let  z  be  the  vertical  direction,  with  the  water  surface  at  z  -  h/2 
and  the  channel  bottom  at  z  -  -h/2  ;  and  let  u' (z)  be  a  z-dependent  pertur¬ 
bation  of  the  depth-averaged  velocity  u  .  By  definition,  u'(z)  must 
satisfy  the  constraint 

h/2 

f  u^ (z) dz  =  0  (1^) 

-K/2 

If  u  +  a'  is  substituted  for  u  when  deriving  the  depth-averaged  momentum 
equation,  then  there  is  a  perturbation  hr'  of  the  depth-integrated  shear 
stress  hr  12  in  the  vertical  plane  (elevation  view) ,  given  by 

h/2 

hr'  =  -p  f  u^Vdz  (15) 

-ff/2 

where  u'  and  v'  are  the  x-  and  y-components  of  u' . 

Equation  14  requires  that  the  depth  integral  be  zero  for  both  u'  and 
v'  ,  but  not  for  their  product  u'v'  .  The  simplest  functional  forms  that 
make  the  integral  of  u'  zero,  without  necessarily  making  the  integral  of 
u'v'  zero,  are  the  expressions 


u'  =  WjZ 

(16) 

V  =  -w^z 

(17) 

where  Wj  and  W2  are  the  x-  and  y-components  of  depth-averaged  vorticity. 

If  the  primary  flow  is  in  the  x-direction,  so  that  v  -  0  but  u  0  , 
then  is  the  streamwlse  vorticity  associated  with  out-of-plane  (secondary) 
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notion  of  fluid  particles  in  the  (primary)  xz-plane.  The  other  vorticity 
component  (i>2  is  perpendicular  to  and  characterizes  the  z-dependence  of 

u' (z)  in  the  xz  plane.  Substitution  of  Equations  16  and  17  into  Equation  15 
produces 

,  _  ph2wiW2  (18) 

■  — n~ 

The  dominant  term  in  the  empirical  relation  used  by  Johannesson  and  Parker 
(1989,  a  and  b)  Indicates  that  a  first  approximation  for  u>2  is 

where  C2  is  a  constant  of  proportionality.  Equations  18  and  19  now  give  way 
to 


<20) 

Assuming  that  Equation  19  holds  for  streamwise  velocity  in  any  direc¬ 
tion,  then  Wi  can  be  equated  with  the  streamwise  vorticity  w,  and  r'  can 
be  eqxiated  with  r,  .  It  is  convenient  to  express  the  latter  as 

r.  =phn|u|^  (21) 

where 


Q 


C2", 

~Tr 


(22) 


With  Cf  given  by  Equation  7,  only  Q  remains  to  be  determined  in  the 
expression  for  r,  . 

Consider  a  vertical  column  of  water  in  a  flow  that  has  (lateral)  radius 
of  curvature  r  in  the  (plan  view)  xy-plane.  Let  the  streamwise  velocity 
\i,(z)  be  given  by 

U,(Z)  «  |ul|l  ♦C2|;crj 

with  the  outward  radial  velocity  Ux(z)  given  by 

Uj(z)  =  w,z 

At  any  vertical  position  z  ,  there  is  a  centrifugal  (outward  radial) 


(23) 


(24) 
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acceleration  u,Vr  due  to  curvature.  This  creates  an  out-of-plane  angular 
acceleration  zxi^^/r  along  the  column,  which  in  turn  produces  the  secondary 
flow. 

To  satisfy  the  no-slip  condition  at  the  very  bottom  of  the  column,  the 
radial  velocity  u^  must  go  from  Wgh/2  to  zero  in  some  small  vertical  dis¬ 
tance  S  .  Assuming  there  is  a  vertical  eddy  viscosity  proportional  to 

u!  S  ,  then  there  must  be  a  shear  stress  at  the  bottom,  roughly  propor¬ 
tional  to  uj  ,  which  opposes  out-of-plane  rotation.  This  stress  is 

negligible  at  distances  greater  than  6  from  the  bottom.  Conservation  of 
angular  momentum  for  the  entire  column  requires  that 

h/2  h/2  2 

J  WaZ^dz  =  |u|  +  p  J  ^zdz  (25) 

-H/2  -H/2 

where  is  a  constant  of  proportionality. 

When  the  integration  is  carried  out  with  Equation  23  substituted  for 
u,  ,  Equation  25  gives  way  to 

^  -  2C,/5r-!|^'  -  (26) 

After  Equations  22  and  26  are  combined,  the  equation  for  streamwise  vorticity 
reduces  to 

dn  ^  c|  ^  \u\^  _  ^  jul  (27) 

Equation  27  asserts  that  dO/dt  depends  upon  velocity,  depth,  bottom 
friction,  radius  of  curvature,  and  two  free  parameters,  and  C2  .  The 

functional  form  comes  from  the  assumed  z-dependence  of  the  streamwise  and  out- 
of-plane  velocities,  and  from  the  imposition  of  centrifugal  and  frictional 
forces  on  a  vertical  column  of  water.  In  this  context,  the  coefficients 
and  C2  are  universal  constants  that  have  to  be  determined  from  experimental 
data. 

Since  Equation  27  is  based  solely  on  qualitative  arguments,  the  grouping 
of  coefficients  is  arbitrary.  The  Important  thing  is  that  the  equation 
contains  two  free  parameters.  Furthermore,  in  addition  to  production  and 
dissipation,  one  also  expects  some  lateral  diffusion  of  n  caused  by  small- 
scale  turbulence.  Based  on  these  considerations,  the  following  is  proposed  as 


12 


the  governing  equation  for  streamwlse  vortlclty: 

+  u.vo  - +  ^v.(^hvn)  (28) 

^  "  rh{l*9hVr2)  '  nr  h  ' 

This  Is  the  secondary  flow  equation  that  Is  used  In  STREMR.  The  coeffi¬ 
cients  and  C2  have  been  replaced  for  convenience  by  A,  and  D,  ,  and 

the  production  term  has  been  reduced  by  an  extra  factor  (1  +  9h^/r^)  In  the 
denominator.  The  latter  modification  leaves  the  production  term  unchanged  as 
long  as  h/r  «  1  .  Otherwise,  the  arguments  supporting  Equation  27  are 
Invalid,  and  production  of  0  Is  automatically  forced  to  zero  when  h/r 
becomes  large. 

With  the  governing  equation  established,  only  the  question  of  initial 
and  boundary  conditions  remains .  The  Initial  value  of  Cl  Is  taken  to  be  zero 
everywhere,  and  Cl  is  fixed  at  zero  for  all  time  on  inlets.  The  normal 
derivative  In  the  diffusion  term  Is  set  to  zero  on  outlets  and  sidewalls. 

Since  the  secondary  flow  must  vanish  on  the  sidewalls,  production  of  Q  Is 
reduced  by  one  half  In  grid  cells  adjacent  to  the  sidewalls.  Otherwise, 

STREMR  discretizes  and  solves  Equation  28  in  the  same  manner  as  Equations  11 
and  12. 
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PART  IV:  RESULTS 


The  secondary  flow  model  (Part  III,  Equation  28)  contains  two  empirical 
coefficients  (A,  and  D,)  whose  values  must  be  established  by  trial-and-error 
comparison  of  predicted  and  measured  velocities.  If  the  SFC  is  genuinely 
reliable,  then  the  values  of  A,  and  D,  for  one  channel  and  one  flow 
condition  should  give,  without  further  adjustment,  satisfactory  predictions 
for  other  channels  and  other  flow  conditions.  It  will  now  be  demonstrated 
that  coefficients  tuned  for  a  single  270-deg*  bendway  can  indeed  be  used  to 
predict  velocity  distributions  for  two  dissimilar  channels  with  multiple  bend- 
ways.  That  does  not  guarantee  that  the  particular  values  used  for  A,  and 
D,  will  be  the  best  for  all  channels  of  interest,  but  it  is  a  crucial  step 
toward  validating  the  SFC.  If  STREMR  were  able  to  reproduce  only  the  bendway 
data  for  which  A,  and  D,  were  specifically  tuned,  then  it  would  Indicate 
that  something  Important  is  still  missing  in  the  secondary  flow  model. 

270-Deg  Bendwav 

The  benchmark  for  establishing  preliminary  values  of  A,  and  D,  is  a 
27()-deg  bendway  with  a  vertical  inner  sidewall  and  a  sloped  outer  bank.  Fig¬ 
ure  1  shows  the  wetted  cross  section,  with  the  STREMR  grid  (387  cells  long  x 
32  cells  wide)  superposed  on  the  plan  view.  All  linear  dimensions  are  given 
in  feet.  For  the  test  under  consideration  (Hicks,  Jin,  and  Steffler  1990), 
the  flow  rate  is  0.83  cfs  and  the  Hanning  coefficient  is  0.010,  which  makes 
Cf  •»  0.0025  in  the  middle  of  the  channel.  The  ratio  of  depth  to  radius  of 
curvature  (h/r)  varies  from  0.016  at  the  toe  of  the  outer  bank,  to  0.021  at 
the  inner  sidewall.  The  outer  bank  was  modeled  in  STREMR  with  four  one-cell- 
wide  stair  steps.  The  inflow  velocity  was  assumed  uniform  except  on  the  outer 
bank,  where  it  was  specified  as  a  linear  function  of  distance  from  the  water's 
edge.  After  considerable  trial  and  error,  the  SFC  coefficients  were  set  at 
A,  -  5  and  D,  -  1/2. 

Figure  2  shows  velocity  vectors,  computed  with  and  without  the  SFC,  at 
stations  along  the  full  length  of  the  channel.  Figure  3  compares  predicted 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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INSIDE  OUTSIDE 
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Predicted  with  SFC 


b.  Predicted  without  SFC 

Figure  2.  Computed  velocity  vectors  for  270-deg  bendway 


a.  Station  1 


Figure  3.  Depth-averaged  velocity  for  the  270-deg  bendway 

(Continued) 
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and  observed  velocities  at  four  individual  stations  (0,  90,  180,  and  270  deg). 
Normalized  velocity  is  defined  to  be  local  depth-averaged  velocity  divided  by 
average  velocity  for  the  channel  cross  section.  Normalized  lateral  position 
is  radial  distance  from  the  inner  sidewall  divided  by  cross-channel  width  of 
the  water  surface.  For  each  station,  a  plan  view  of  the  channel  is  provided 
with  an  arrow  that  indicates  the  location  of  the  data  station  along  the  inner 
sidewall  of  the  bendway.  The  velocity  profiles  are  shown  as  they  would  be 
seen  by  an  observer  looking  in  the  direction  of  the  arrow. 

The  SFC/STREMR  velocity  profiles  reproduce  the  observed  profiles  well 
except  on  the  outer  bank,  where  the  computed  bottom  resistance  is  consistently 
too  low.  The  underprediction  also  occurs  in  straight  channels,  and  it 
indicates  the  need  for  some  modification  of  Equation  7  in  response  to  bottom 
slope.  Finer  discretization  (more  grid  lines  concentrated  on  the  bank)  does 
not  eliminate  the  problem.  It  remains  to  be  seen  whether  this  deficiency  can 
be  cured  without  using  a  fully  3-D  numerical  model.  Otherwise,  the  SFC  (with 
its  tuned  coefficients)  does  a  good  job  of  making  STREMR  predictions  match  the 
test  data. 

Without  the  SFC,  the  hipest  velocity  remains  near  the  inside,  and  the 
computed  profiles  bear  little  resemblance  to  the  observed  profiles  downstream 
of  the  bend  entrance  (Station  1).  In  contrast,  the  SFC  causes  high  velocity 
to  migrate  to  the  outside  because  of  the  production  term  in  Equation  28,  which 
becoBKS  nonzero  whenever  h/r  is  nonzero.  The  production  term  returns  to 
zero  in  straight  sections,  leaving  only  the  dissipation  term,  which  gradually 
kills  the  secondary  flow.  If  there  is  a  reversal  in  curvature,  there  will  be 
a  migration  of  high  velocity  toward  the  opposite  side  of  the  channel.  The 
coefficients  A,  and  D,  determine  the  precise  rates  of  migration,  dissipa¬ 
tion,  and  reversal. 


Channel  Bend  Facility 

Maynord*  has  made  detailed  velocity  measurements  in  an  S-shaped  flume 
called  the  Channel  Bend  Facility  (CBF) .  The  CBF  entails  two  bends  with  a 
reversal  in  curvature  and  a  trapezoidal  cross  section,  and  it  represents  a 


*  Unpublished  test  data  provided  by  S.  T.  Maynord,  Research  Hydraulic 
Engineer,  August  1987,  US  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  NS. 
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considerable  departure  from  the  geometry  of  the  270-deg  bendway.  Figure  4 
shows  the  wetted  cross  section,  with  the  STREMR  grid  (121  cells  long  x 
46  cells  wide)  superposed  on  the  plan  view.  All  linear  dimensions  are  given 
in  feet.  For  the  test  under  consideration,*  the  flow  rate  is  6.75  cfs  and 
the  Manning  coefficient  is  0.02,  which  makes  -  0.0075  in  the  middle  of 
the  channel.  The  ratio  of  depth  to  radius  of  curvature  (h/r)  varies  in  the 
bends  from  0.018  at  the  toe  of  the  outer  bank,  to  0.025  at  the  toe  of  the 
inner  bank.  The  banks  were  discretized  in  STREMR  with  five  one-cell-wide 
stair  steps  on  each  side  of  the  channel.  The  inflow  velocity  was  assumed 
tiniform  except  on  the  banks,  where  it  was  specified  as  a  linear  function  of 
distance  from  the  water's  edge.  The  SFC  coefficients  were  set  at  A,  -  5  and 
D,  -  1/2. 

Figure  5  shows  velocity  vectors,  computed  with  and  without  the  SFC,  at 
stations  along  the  full  length  of  the  CBF.  Figure  6  compares  predicted  and 
observed  velocities  at  five  individual  stations,  whose  locations  are  indicated 
by  the  arrow  on  each  included  plan  view.  As  before,  the  velocity  profiles  are 
shown  as  they  would  be  seen  by  an  observer  looking  in  the  direction  of  the 
arrow.  Normalized  velocity  is  defined  to  be  local  depth-averaged  velocity 
divided  by  average  velocity  for  the  channel  cross  section.  Normalized  lateral 
position  is  radial  distance  from  the  inner  water  surface  (indicated  by  the  tip 
of  the  arrow)  divided  by  cross-channel  width  of  the  water  surface. 

The  agreement  between  the  SFC/STREMR  predictions  and  the  observed 
velocity  profiles  is  at  least  as  good  for  the  CBF  as  it  was  for  the  270-deg 
bendway.  Here  the  average  value  of  h/r  is  only  slightly  greater,  but  the 
value  of  is  three  times  greater  here  than  that  for  the  270-deg  bend. 

This  lends  support  to  the  proposed  influence  of  C^  in  the  SFC. 

Riprap  Test  Facility 

The  Riprap  Test  Facility  (RTF)  is  a  trapezoidal  channel  with  four  bends 
and  two  reversals  in  curvature.  Maynord*  has  conducted  a  test  in  the  RTF  with 
a  flow  rate  of  49.5  cfs  and  a  Manning  coefficient  of  0.026.  Figure  7  shows 
the  wetted  cross  section,  with  the  STREMR  grid  (391  cells  long  x  36  cells 


*  Unptdalished  test  data  provided  by  S.  T.  Maynord,  Research  Hydraulic 
Engineer,  August  1987,  US  Aimy  Engineer  Waterways  Experiment  Station, 
Vicksburg,  MS. 


20 


CHANNEL  CROSS  SECTION 


Grid  plan  view  and  channel  cross  section  for  CBF 


Predicted  with  SFC 


b.  Predicted  without  SFC 
Figure  5.  Computed  velocity  vectors  for  CBF 
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Figure  6.  (Sheet  2  of  3) 
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e .  Station  5 

Figure  6.  (Sheet  3  of  3) 
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17.2 


Figure  7.  Grid  plan  view  and  channel  cross  section  for  RTF 


wide)  superposed  on  the  plan  view.  All  linear  dimensions  are  given  In  feet. 
The  value  of  Cf  is  0.0088  In  the  middle  of  the  channel,  and  h/r  varies  in 
the  bends  from  0.024  at  the  toe  of  the  outer  bank  to  0.031  at  the  toe  of  the 
inner  bank.  The  average  value  of  h/r  is  about  50  percent  greater,  and  the 
value  of  Cf  is  3.5  times  greater  here  than  that  for  the  270-degree  bend. 

The  banks  were  discretized  in  STREMR  with  six  one-cell-wide  stair  steps  on 
each  side  of  the  channel.  The  inflow  velocity  was  assumed  uniform  except  on 
the  banks,  where  it  was  specified  as  a  linear  function  of  distance  from  the 
water's  edge.  The  SFC  coefficients  were  set  at  A,  -  5  and  D,  -  1/2. 

Figure  8  shows  velocity  vectors,  computed  with  and  without  the  SFC,  at 
stations  along  the  full  length  of  the  RTF.  Figure  9  compares  predicted  and 
observed  velocities  at  13  individual  stations,  whose  locations  are  indicated 
by  the  arrow  on  each  included  plan  view.  As  before,  the  velocity  profiles  are 
shown  as  they  would  be  seen  by  an  observ^er  looking  in  the  direction  of  the 
arrow.  Normalized  velocity  is  defined  to  be  local  depth-averaged  velocity 
divided  by  average  velocity  for  the  channel  cross  section.  Normalized  lateral 
position  is  radial  distance  from  the  inner  water  surface  (indicated  by  the  tip 
of  the  arrow)  divided  by  cross-channel  width  of  the  water  surface. 

Up  to  the  third  bend,  the  agreement  between  the  SFC/STREMR  predictions 
and  the  observed  velocity  profiles  is  about  as  good  for  the  RTF  as  it  was  for 
the  CBF  and  the  270-deg  bendway.  There  is  some  deterioration  of  accuracy  in 
the  third  bend  (Stations  10  through  12) ,  even  though  the  predictions  still 
follow  the  observed  trends  for  velocity  migration.  Accuracy  has  begun  to 
recover  at  the  entrance  to  the  fourth  bend  (Station  13) . 

Hypothetical  Effect  of  Bend  Angle 

Since  the  SFC  (with  A,  -  5  and  D,  -  1/2)  yields  acceptable  predic¬ 
tions  for  three  rather  different  channels,  one  expects  that  it  should  be  no 
less  applicable  for  minor  variations  on  these  configurations.  With  this  in 
mind.  STREMR  was  used  to  calculate  the  hypothetical  effect  of  bend  angle  for  a 
single-bend  channel  with  the  same  cross  section,  radius  of  curvature,  and 
friction  coefficient  as  the  CBF  (Figure  4).  Specifically,  unverified  predic¬ 
tions  were  made  for  bend  angles  of  30,  60,  90,  and  120  deg,  with  20-ft 
straight  sections  at  the  entrance  and  exit  of  the  bend.  The  grid  spacing  was 
the  same  as  that  for  the  CBF. 
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a.  Predicted  with  SFC 


Figure  8.  Computed  velocity  vectors  for  RTF 


Figure  9.  Depth-averaged  velocity  for  the  RTF  (Sheet  1  of  7) 


w 


29 


d.  Station  4 
Figure  9.  (Sheet  2  of  7) 
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f.  Station  6 
Figure  9.  (Sheet  3  of  7) 
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Morm*hx«il  Velocity 


g.  Station  7 


Normoitsod  Velocity 

h.  Station  8 
Figure  9.  (Sheet  4  of  7) 
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j .  Station  10 
Figure  9.  (Sheet  5  of  7) 
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NormatiMd  VcJoctljr 


k.  Station  11 


1.  Station  12 
Figure  9.  (Sheet  6  of  7) 
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m.  Station  13 
Figure  9.  (Sheet  7  of  7) 
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Figures  10  through  13  show  velocity  vectors  computed  with  and  without 
the  SFC  for  the  four  different  bend  angles .  As  in  previous  cases ,  the  SFC 
moves  the  high  velocities  gradually  toward  the  outside  of  the  bend,  where  they 
remain  for  some  distance  downstream.  Its  effect  is  noticeable  even  for  the 
30-deg  bend.  In  contrast,  omission  of  the  SFC  keeps  the  high  velocity  near 
the  insides  of  the  bends. 
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b.  Predicted  without  SFC 


Computed  velocity  vectors  for  hypothetical  30-deg  bend 


b.  Predicted  without  SFC 

Figure  11.  Computed  velocity  vectors  for  hypothetical  60-deg  bend 
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b.  Predicted  without  SFC 

Figure  12.  Computed  velocity  vectors  for  hypothetical  90-deg  bend 
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b.  Predicted  without  SFC 


Figure  13.  Computed  velocity  vectors  for  hypothetical  120-deg  bend 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


An  empirical  governing  equation  has  been  proposed  for  the  secondary  flow 
that  gives  rise  to  migration  of  high  velocity  toward  the  outside  of  channel 
bends.  The  secondary  flow  stems  from  the  interaction  between  lateral  curva¬ 
ture  and  vertically  nonuniform  velocity,  and  it  creates  an  unbalanced  force 
that  alters  the  primary  flow.  The  depth-averaged  secondary  flow  correction 
(SFC)  Includes  two  coefficients  (A,  and  D,)  whose  values  have  been  tuned  for 
agreement  between  STREMR  numerical  model  predictions  and  experimental  data  for 
a  270-deg  bendway.  Using  the  same  values  for  A^  and  D,  ,  the  SFC  yields 
STREMR  predictions  with  comparable  accuracy  for  the  CBF  (two  bends)  and  the 
RTF  (four  bends),  both  of  which  have  depth,  bottom  friction,  and  radius  of 
curvature  different  from  the  270-deg  benchmark. 

Results  obtained  with  the  SFC  are  encouraging,  since  it  was  not  clear  in 
advance  that  the  same  values  of  A,  and  D,  would  render  acceptable  predic¬ 
tions  for  more  than  a  single  channel  configuration.  The  form  of  the  secondary 
flow  equation  ensures  proper  qualitative  behavior,  but  not  necessarily  quanti¬ 
tative  accuracy.  Comparisons  of  predicted  and  measured  velocities  have  demon¬ 
strated,  however,  that  coefficients  tuned  for  one  channel  can  also  be  used  for 
others.  It  appears  that  the  SFC  may  be  used  (with  A,  -  5  and  D^  -  1/2) 
for  ratios  of  depth  h  to  radius  of  curvature  r  in  the  range  0.0  <  h/r  < 
0.04  ,  and  for  friction  coefficients  in  the  range  0.002  <  Cf  <  0.01  .  It 
remains  to  be  seen  how  well  the  SFC  works  outside  these  bounds  without  chang¬ 
ing  the  values  of  A,  and  D,  .  To  answer  the  latter  question,  more  tests 
will  be  needed  for  other  channels  and  flow  conditions. 
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APPENDIX  A:  NOTATION 


A. 

Ci,C2,C^ 

C1.C2 

Cf 

D. 

h 

k 

n 

n 

p 

r 

£ 

t 

I 

u 

a 

u' 

U. 

vi/x 

U' 

V 

V' 

X 

z 


Empirical  production  coefficient  for  secondary  flow 
Empirical  coefficients  in  turbulence  model 
Constants  of  proportionality 
Friction  factor  given  by  Manning's  equation 
Empirical  decay  coefficient  for  secondary  flow 
Depth 

Turbulence  energy 
Manning's  coefficient 
Unit  vector  normal  to  u 
Pressure 

Lateral  radius  of  curvature 
Force  arising  from  secondary  flow 
Time 

Viscous  force  arising  from  the  depth-averaged  stress  tensor 

x-component  of  u 

Depth— averaged  vector  velocity 

x-component  of  jil 

Radial  velocity 

Outward  radial  velocity 

Streamwise  velocity 

Centrifugal  (outward  radial)  acceleration 
z-dependent  perturbation  of  u 
y— component  of  u 
y-component  of  u' 

Resistance  force  (per  unit  mass)  due  to  bottom  friction 
Vertical  position 


A1 


zu2/r 

e 

V 

P 

r. 


«2 


Q 

V 


Out-of-plane  angular  acceleration 
Turbulence  dissipation  rate 
Eddy  viscosity 
Density 

Depth-averaged  shear  stress 

X— component  of  depth-averaged  vorticity 

y— component  of  depth-averaged  vorticity 

Streamwise  component  of  depth-averaged  vorticity 

C  w 
_2_» 

12 

Gradient  operator 


A2 
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